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ABSTRACT 

We present a measurement of the absolute surface brightness of the zodiacal light (3900-5I00A) 
toward a fixed extragalactic target at high ecliptic latitude based on moderate resolution (^1.3A per 
pixel) spectrophotometry obtained with the du Pont 2.5m telescope at Las Campanas Observatory in 
Chile. This measurement and contemporaneous Hubble Space Telescope data from WFPC2 and EOS 
comprise a coordinated program to measure the mean flux of the diffuse extragalactic background light 
(EBL). The zodiacal light at optical wavelengths results from scattering by interplanetary dust, so that 
the zodiacal light flux toward any extragalactic target varies seasonally with the position of the Earth. 
This measurement of zodiacal light is therefore relevant to the specific observations (date and target 
field) under discussion. To obtain this result, we have developed a technique that uses the strength of 
the zodiacal Fraunhofer lines to identify the absolute flux of the zodiacal light in the multiple-component 
night sky spectrum. Statistical uncertainties in the result are 0.6% (Icr). However, the dominant source 
of uncertainty is systematic errors, which we estimate to be 1.1% (Icr). We discuss the contributions 
included in this estimate explicitly. The systematic errors in this result contribute 25% in quadrature to 
the final error in our coordinated EBL measurement, which is presented in the first paper of this series. 

Subject headings: Diffuse radiation — cosmology: observations — techniques: spectroscopic — 
interplanetary medium 

1. INTRODUCTION 

This is the second in a series of three papers in which 
we present a measurement of the mean flux of the opti¬ 
cal extragalactic background light and the cosmological 
implications of that result (see Bernstein, Freedman, & 

Madore 2002a & 2002c). The extragalactic background 
light (EBL) is the spatially averaged surface brightness of 
all extragalactic sources, resolved and unresolved. As such, 
the absolute flux of the EBL is a powerful and fundamen¬ 
tal cosmological constant which can significantly constrain 
galaxy formation and evolution scenarios. Like all diffuse 
backgrounds, however, the optical EBL is very difficult to 
isolate from foreground sources, which are two orders of 
magnitude brighter. At high Galactic and ecliptic lati¬ 
tudes (> 30°), the sky flux observed from the ground is 
dominated by terrestrial airglow and zodiacal light (ZL), 
each with a surface brightness of ~ 23 AB mag arcsec“^. 

The Hubble Space Telescope (HST), which orbits at an 
altitude of 600 km, avoids atmospheric emission, but the 
total sky flux is still dominated by ZL. An accurate mea¬ 
surement of the ZL is therefore crucial to a successful de¬ 
tection of the diffuse EBL. 

Our measurement of the EBL involves simultaneous 
HST and ground-based observations. From HST we mea¬ 
sure the total flux of the night sky, including ZL. Using 
spectrophotometry over the range 3860-5150A (1.25A per 
pixel) taken with the Boiler & Ghivens Spectrograph on 
the duPont 2.5m Telescope at Las Gampanas Observatory 
in Chile, we measure the absolute flux of the ZL contribut¬ 


ing to the HST observations, which we can then subtract. 
In Bernstein, Freedman, & Madore (2002a, henceforth Pa¬ 
per I), we present the full details of the coordinated pro¬ 
gram to measure the EBL. In this paper, we present the 
ground-based measurement of the absolute flux of the ZL. 
As calibration of these data is crucial to the scientific goals, 
the data acquisition, reduction, and flux calibration are 
discussed here in detail. 

Background regarding the nature of the ZL is given in 
§^. The observations, data reduction, and flux calibra¬ 
tion are discussed in §|[ In §^, we briefly described the 
complications which arise due to atmospheric scattering, 
which redirects off-axis flux into and on-axis flux out of 
the line of sight. Detailed calculations of the atmospheric 
scattering relevant to precisely our observing situation (de¬ 
fined by the observatory location and positions of the Sun, 
Galaxy, and target relative to eachother and the horizon) 
are relegated to the Appendix, and summarized in § ^ In 
§^, we describe the technique used to measure the zodiacal 
light flux in reduced spectra. The results are summarized 
in §1 

2. BACKGROUND 

Zodiacal light (ZL) is sunlight scattered off of dust 
grains in the solar system. Dust column densities are 
largest toward the ecliptic plane, causing the highest ZL 
intensities there. The scattering geometry is illustrated in 
Figure]^ for a high latitude field viewed near the anti-solar 
direction. The figure shows the scattering angle and line of 
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sight as defined by geocentric ecliptic longitude (A —A©, in 
which A is ecliptic longitude) and ecliptic latitude (/3). At 
ecliptic latitudes greater than 30 degrees and large scatter¬ 
ing angles, the zodiacal light can be as faint as ~ 23 mag 
arcsec”^ at SOOOA (^1x10“^ ergs s“^ cm“^ sr“^ A“^) 
with a solar-type spectrum. In the ecliptic plane and near 
the Sun, the ZL can be as much as 20 times brighter than 
at high ecliptic latitudes and significantly reddened rela¬ 
tive to the incident solar spectrum. The positional de¬ 
pendence of the ZL flux and color is a function of the 
density and composition of the scattering particles and of 
scattering geometry such that the zodiacal light is redder 
and brighter at smaller elongation angles (see Leinert et al. 
1998). In general, the solar spectrum is preserved in the 
spectrum of the ZL with less than 30% deviation in the 
broad-band spectral shape from the UV to the near-IR 
(0.2 — 2/im). 

Such weakly wavelength-dependent scattering can be 
characterized by simple Mie scattering theory, which gen¬ 
erally describes the interaction of photons with solid parti¬ 
cles larger than the wavelength of the incident light (Roser 
& Staude 1978). The interplanetary dust (IPD) cloud at 
the orientations of interest to us is known to be composed 
predominantly of particles larger than 10/rm with mod¬ 
erate surface roughness and layered composition. These 
conclusions are based on IR observations of the thermal 
properties of the IPD, dynamical arguments, and labora¬ 
tory work on dust captured in the upper atmosphere and 
on the moon (see Reach et al. 1996; Berriman et al. 1994, 
Dermott etal. 1996; Brownlee 1978; Fechtig etal. 1974; 
Leinert etal. 1998 and references therein). In the absence 
of free electrons, the cross-sections for non-linear scatter¬ 
ing, such as Raman scattering or two-photon processes, 
are too low to be a significant effect for dust particles with 
large dielectric constants. 

Indeed, Mie scattering models for rough particles with 
the size and composition of the IPD successfully describe 
the weakly wavelength-dependent scattering characteris¬ 
tic of the ZL: the incident spectrum is slightly reddened 
over very broad band-passes, while narrow-band spec¬ 
tral features are preserved to high precision (Weiss-Wrana 
1983, Schiffer 1985). Both these narrow- and broad-band 
characteristics have been empirically demonstrated, as the 
ZL spectrum is seen to be roughly 5-10% redder per lOOOA 
over the range from 3000-8000A than the Sun (see Leinert 
etal. 1998 and references therein), and Beggs etal. (1964) 
have shown that the Fraunhofer lines in the ZL show no 
measurable deviation from their solar equivalent widths to 
the accuracy of their calibration (2%). The broad-band 
reddening is characterized as the color of the ZL relative 
to the solar spectrum as a function of wavelength, 


C(A,Ao) 


Izl{X)IIq{\) 

Izl{Xq)/Iq{\o) 


( 1 ) 


While Mie scattering explains the spectrum and inten¬ 
sity of the ZL in general, the exact flux of the ZL cannot 
be modeled to the accuracy we require here, simply due 
to uncertainties in the exact composition and column den¬ 
sity of the IPD. Empirically, the surface brightness of the 
ZL is known to roughly 10% accuracy as a function of 
scattering geometry and ecliptic latitude alone. However, 
uncertainty in flux calibrations of those surface bright¬ 
ness measurements and spatial variability in the IPD cloud 



Fig. 1. — The geometry of zodiacal light scattering (adapted from 
Matsuura, Matsumoto, &; Matsuhara 1995). The viewing line of sight 
is defined by ecliptic latitude (/3), and geocentric ecliptic longitude, 
(A — Aq), or alternatively, by f) and the elongation angle, e, which is 
defined as cos e = cos(A — Aq) cos /3. The scattering angle is 9. 


(such as cometary trails) preclude more accurate predic¬ 
tions (Levasseur-Regourd & Dumont 1980; Richter etal. 
1982; Leinert etal. 1998 and references therein). In addi¬ 
tion, most published measurements of the ZL flux include 
the background EBL in the reported surface brightness of 
the ZL. In order to detect the EBL, we must explicitly mea¬ 
sure the zodiacal contribution to the sky surface brightness 
along a specific line of sight at a particular wavelength. 

In order to uniquely define the ZL spectrum, both the 
mean flux in a narrow-bandpass and the broad-band color 
over the wavelength range of interest must be measured for 
the line of sight in question. We can do so by making use 
of the fact that the equivalent widths of solar Fraunhofer 
lines, reproduced in the ZL spectrum, are known to very 
high accuracy: we can then determine the continuum level 
(mean flux) of the ZL at a given wavelength by measuring 
the apparent equivalent width of the Fraunhofer in spec¬ 
trum of the night sky. For example, if the equivalent width 
of the Fraunhofer lines in the night sky around 4500A are 
only one half the strength of the same absorption features 
in the solar spectrum, then the ZL contributes only 50% of 
the night sky flux at that wavelength. A practical compli¬ 
cation of this approach is, of course, identifying the con¬ 
tinuum level for accurate measurement of the equivalent 
widths. Because the ZL and atmospheric emission (air- 
glow) contribute almost equally to the night sky flux in 
the visible wavelength range at the ecliptic orientation of 
our observations, this becomes prohibitively difficult where 
airglow lines become strong (above 5500A). The problem is 
further complicated by the fact that the airglow spectrum 
is continually changing, both in the strength of particular 
lines and in its mean, wide-band flux. In order to ac¬ 
complish this measurement, we have therefore developed 
a new technique for deconvolving one known signal from 
a variable, multicomponent spectrum. The details of this 
technique are discussed in §|^. 

3. OBSERVATIONS AND DATA REDUCTION 

We obtained long-slit spectra of the night sky using 
the Boiler & Chivens Spectrograph on the 2.5m duPont 
Telescope at Las Campanas Observatory (LCO) on the 
nights of 1995 November 26-29. The first night of the run 
was lost to clouds. On the third night, a thermal short 
in the dewar caused poor temperature regulation of the 
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Fig. 2.— The slit position during program observations. The 
image is a 4 x 4 arcmin^, r-band exposure taken the 2.5m duPont 
telescope at Las Campanas Observatory. WFPC2/HST observations 
were taken in parallel with the spectroscopy presented here for the 
purpose of measuring the EBL. For comparison, the WFPC2/HST 
field of view is also contained within this image, located towards the 
lower left quadrant (see Figure 7, Paper I). 


CCD which re sult ed in low and variable charge transfer 
efficiency (see §3.1). The data from these two nights were 
therefore dropped from further analysis. The dewar was 
repaired for the rem ainin g nights of the run. As discussed 
further in §3.3 and (:3.£, both detector performance and 
photometric conditions were stable to better than 1% on 
the two useful nights of the run. 

As described in Paper I, these spectra were taken within 
the field of view of our HST/WFPC2 observations, which 
executed in consecutive orbits on 1995 November 27-28. 
Because the HST and LCO observations executed simulta¬ 
neously and along exactly the same pointing (see Paper I), 
the ZL contribution to both data sets is identical. The slit 
position is shown overlayed on a r-band image of the field 
in Figure |[ The exact coordinates of the spectroscopic 
observations were selected from ground-based imaging to 
avoid objects brighter than r = 26 mag arcsec“^. This 
was done for convenience in the data reduction; extra- 
galactic sources need not be rigorously avoided, as spec¬ 
tral lines from extragalactic sources fainter than r ~ 24 
AB mag will be significantly redshifted, on average, and 
should not align with solar features. In addition, such ob¬ 
jects will have low enough surface number density (< 0.2 
arcsec”^) that they will not significantly impact the aver¬ 
age extracted spectrum. 

We used a 600 1/mm grating to obtain spectra over the 
wavelength range 3860-5150A with a dispersion of ~ 1.3A 
per pixel. The wavelength range was chosen to include a 
maximum number of strong Fraunhofer lines while avoid¬ 
ing strong airglow features. The strong Mgi Fraunhofer 
lines near ~ 5170A were identified in an earlier run to 
be affected by rapidly variable airglow features and were 
therefore excluded from our 1995 spectral coverage. Ca H 
& K were included in our observations at the blue end. 


A slit-width of ~ 1.5arcsec (see ( j3.6.3D produced roughly 
2.6A resolution in the program observations. Even though 
the ZL has a surface brightness of fiy ~23.2 mag arcsec“^, 
we obtained a signal-to-noise ratio of ~ 40 per spectral 
resolution element from a single, 30 minute exposure by 
integrating over the total slit surface area (^300arcsec^). 
To minimize read-noise, we binned the data on-chip by 
four pixels in the spatial direction and averaged over the 
full 3.4 arcmin spatial extent of the slit in the data reduc¬ 
tion. 

Spectrophotometric standards were observed through a 
slit 10.8arcsec wide. The detector we used is a thinned 
Tektronics/SITe CCD with 1024 x 1024 x 24/im pixels, 
1.15e“/DN gain, and 8.6e“ read-noise in the 4x1 binned 
configuration. The quantum efficiency of the CCD is near 
50% over the spectral range of these observations; how¬ 
ever, the spectrograph throughput drops by a factor of 
two between 5000A and 3800A, as can be seen from the 
sensitivity curve plotted in Figure The upper panel of 
Figure || shows the wavelength calibrated spectrum ob¬ 
tained from one spatial resolution element (1 column) of 
one of the 4x1 binned exposures; the lower panel shows 
the averaged spectrum from the full image (87 columns). 
Above 4100A, the count rate from sky is more than twice 
the dark rate. The maximum error is 10% per resolution 
element at the blue end of the spectra, simply due to the 
low count-rate at bluer wavelengths. Above 4100A, the 
error per resolution element is roughly 1%. 

Final calibration errors are summarized in Table |^. The 
data reduction steps are described below in the order 
which they were performed. 

3.1. Detector Linearity 

Because our observations of the ZL have total count lev¬ 
els in the range 20-50 DN pixel” and the standard star 
observations have close to 5000 DN pixel” it is crucial 
to verify that the CCD is linear over this broad range. In 
these data, 16 rows are read off beyond the physical ex¬ 
tent of the chip, averaged together, and recorded as a bias 
row, which can be used as an accurate diagnostic of the 
charge transfer efficiency (CTE) of the CCD. Because the 
slit only illuminates the central third of the chip, any resid¬ 
ual charge which is not passed through the parallel gates 
(due to low CTE) will appear as a jump in the charge 
level of the bias row at the boundary between the exposed 
and unexposed regions of the chip. It is evident from this 
diagnostic that the temperature regulation of the chip be¬ 
came erratic during the third night of the run, causing an 
increase in the spurious charge and causing the CTE to 
drop to an unacceptable level (~ 99.995% per transfer, or 
95.0% over 1024 rows). The data from this night were ex¬ 
cluded from the analysis because charge shared between 
rows due to incomplete charge transfer will affect the ap¬ 
parent strength of spectral features. On the nights during 
which the temperature of the CCD remained stable, the 
mean level in the bias row was not detectably higher in 
columns illuminated by the slit than in those which were 
not illuminated. This was true even for dome flats images, 
in which the mean level of the illuminated columns was 
~7000DN pixel”^. In addition, our final results are based 
on data imaged in rows 1-700 of the CCD. The data which 
contribute to our final results were obtained with a min- 
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Fig. 3.— Signal-to-noise of program data. The upper plot shows the wavelength-calibrated spectrum in electrons per pixel obtained from 
one spatial resolution element (column) of an 1800 second exposure of the night sky. The lower plot shows the spectrum produced by averaging 
87 columns. The corresponding error spectra, both read-noise dominated, are shown with a dotted line at the bottom of each plot. 


imum CTE of 99.9995% per transfer, or 99.7% over 700 
rows. 

Another possible cause of non-linearity at low count lev¬ 
els is deferred charge. We performed the standard tests for 
deferred charge as described in Gilliland (1992) and find 
that only 0.3% of the pixels showed deviations from linear¬ 
ity greater than four times the read-noise. Pixels exhibit¬ 
ing non-linearity were flagged in all images and excluded 
from analysis. We also performed a standard linearity test 
by taking dome flat-held exposures with integration times 
between 0.5 and 200 seconds and looking for variation in 
the detected count-rate between 20-20,000 DN pixel“^. 
The influence of lamp instability was minimized by taking 
several series of exposures and averaging the results. The 
detector response was linear to the limits of the sampled 
range which more than brackets the signal level of the stan¬ 
dard star observations (peaking at roughly 5,000 DN) and 
the program observations (20-50DN on average). Non¬ 
repeating deviations of less than 1% were attributed to 
the instability of the lamp. 

3.2. Bias Subtraction 

As discussed in above, we use the strength of the so¬ 
lar Fraunhofer lines in the ZL spectrum to measure the 


ZL flux. This is a differential measurement in the dis¬ 
persion direction; uniform, additive offsets due to bias or 
dark current will not affect our results beyond the small 
effect on flux calibration. Structure in either bias or dark 
current, however, will increase random errors in the re¬ 
sults. For example, sharp features in the dispersion direc¬ 
tion will increase the rms errors in the average strength of 
the Fraunhofer absorption lines. Also, because we extract 
a single, one dimensional spectrum from every two dimen¬ 
sional image by averaging over the full spatial extent of 
the spectrum, any fluctuations in bias level in the spatial 
direction add random errors to the averaged flux found at 
any wavelength. 

To remove spatial variations in both directions, the bias 
correction was done in three steps. Variations in the 
dispersion direction were subtracted using the 150 col¬ 
umn overscan region. The overscan was fitted with a 
Savitsky-Golay routine (Press etal. 1992), which follows 
rapid jumps in the mean level, and the fit was then sub¬ 
tracted from each column. The mean bias level was then 
removed by subtracting a single mean bias value, which 
is the average of the bias row (row 1024) in each frame. 
Finally, bias variations over the chip were found to be very 
stable after the overscan column and mean bias level were 
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Table 1 

Error Budget for Zodiacal Light Flux 


Step 


Bias level removal (§3.21 _ 

Dark current removal (§ p.3| ) _ 

Pixel-to-pixel flat -fie dine (§BJ) 

Slit illumination (13.4)^ _ 

Point source flux cai [Sp 
Aperture cor rectiqn (§3.6.2|) 

Solid angle (§|3.6.3|) _, 

Lcat(ZL) estimate (§A^ 
Lcat(ISL) estimate (§ST3) 
rms scatter in ZL correlation (§| 


Combined (Ict)® 


Statistical Uncertainty Systematic Uncertainty 


< 0.01% 

< 0.01% 

< 0.01% 

< 0.1% 


0.6% 


0.2% 


0.6% 


1.2% 


0.5% 


0.6%^ 

0 . 6 % 1 . 1 % 


^Statistical errors have been combined in quadrature to obtain a cumulative, one-sigma error. 
Systematic errors have been combined assuming a flat probability distribution for each contribut¬ 
ing source of error. The resulting systematic error is roughly Gaussian distributed, and the quoted 
value is the 68% confidence interval. For a detailed discussion see Paper I. 


subtracted from every image. These variations were thus 
removed by subtracting a “superbias” image, which is the 
average of 250, overscan- and bias-subtracted bias frames. 
To verify the accuracy of this procedure, we test-reduced 
50 bias frames which had not been included in the su¬ 
perbias. After bias subtraction, these test frames had an 
average value of 0.005 DN with an rms scatter of 0.03 DN 
and no residual systematic structure. 

3.3. Dark Current Subtraction 

Dark frames were taken immediately before and after 
observations on each night, and 20 darks were taken at 
the start and end of the run, which were combined to 
make a “superdark,” the mean level of which drifted by 
~ 1.0 DN over the extent of the frame. Unfortunately, the 
20-frame superdark is read-noise dominated and cannot 
provide a pixel-to-pixel correction. The superdark was 
therefore smoothed using a sliding 3x3 boxcar median 
filter, in order to avoid adding noise and to allow the re¬ 
moval of the mean dark level. The frames taken during the 
run were used to test the accuracy of the superdark: after 
bias subtraction (as described above) and dark subtrac¬ 
tion, the test-reduced darks taken at the beginning and 
end of the observations on 1995 November 27 and 29 had 
a mean level of ±0.25 DN with no coherent pattern. 

3.4. Flat Fielding and Illumination Correction 

We used a 1.5 arcsec slit for the program observations 
in order to preserve resolution. However, we used a 10.8 
arcsec slit for the standard star observations in order to 
collect as much light as possible. Different sets of flat- 
field and illumination corrections were therefore required 
for two reasons. First, microscopic roughness on the edges 
of the slit jaws caused shadowing which changed as a func¬ 
tion of slit-width. Second, the slit-jaw mechanism on this 
spectrograph is such that the jaws are parallel for separa¬ 
tions up to ~ 5.4 arcsec (500/rm) but are not parallel when 
the jaws are set to 10.8 arcsec in the center. Variation in 
slit-width is almost 10% from end-to-end when the width 
at the center is 10.8 arcsec. To compensate for the vari¬ 


able slit-width for standards, the illumination corrections 
for both slit-widths were normalized to the spatial cen¬ 
ter of the slit. Standard stars were all observed within 
two pixels of the central column used for normalization, 
which places them within 99.98% of the nominal value for 
the wide-slit illumination correction. Flat-held and slit- 
illumination corrections were created from dome and twi¬ 
light sky flats, respectively, using the tasks RESPONSE 
and ILLUMINATION in the IRAF SPECRED package. 

3.5. Wavelength Calibration 

Wavelength solutions were based on He-Ar comparison 
spectra taken before and after each 1800 second program 
exposure and immediately after each standard star ob¬ 
servation. At the beginning of the run, great care was 
taken to align the dispersion axis with the pixel rows: the 
rms variation in the centroid position of arc lines over 
the full 3.4 aremin spatial extent of the images is typi¬ 
cally 0.05 - 0.12 pixels. This is true over the full wave¬ 
length range. Consequently, it was not necessary to rec¬ 
tify the two-dimensional program spectra of the zodiacal 
light (night sky). Each program spectrum was simply aver¬ 
aged along rows (with 5cr cosmic ray rejection) to obtain 
a one-dimensional spectrum, and a wavelength solution 
based on the spatial center of the image was applied af¬ 
terward. Standard stars were first extracted to produce 
one-dimensional spectra and then wavelength calibrated 
in the usual way. 

We identified 21-25 features in each comparison spec¬ 
trum visually and fitted a third order (four term) Legen¬ 
dre polynomial to the pixel-wavelength solution to obtain 
a dispersion curve. This was done using the IDENTIFY 
task in IRAF. The rms residuals in the wavelength solu¬ 
tion were 0.016-0.04A. Shifts in the wavelength solution 
between program observations were less than ±1.5 pixels 
from the start to the end of the night. Linear dispersion 
solutions were applied using the task DISPCOR. 

3.6. Flux Calibration 
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Fig. 4.— The extinction function derived from observations of 
standard stars taken on 1995 November 27 and 29. Extinction so¬ 
lutions for the two nights individually were found to be identical to 
within the one-sigma statistical errors (< 0.01 mag). The function 
shown is a sixth-order Legendre polynomial fit to extinction coeffi¬ 
cients obtained in 50A wide bins. 


Fig. 5.— Two sensitivity functions obtained from observations of 
tertiary standards on different nights. The solutions differ by 1-2%, 
well within the fluctuations expected to arise from small temperature 
variations of the CCD. 


Three independent components of the surface brightness 
calibration will affect its final accuracy: point source cali¬ 
bration (which includes the sensitivity and extinction cor¬ 
rections); aperture correction (which compensates for the 
loss of light from point source observations which does not 
occur in observation of a uniform aperture-filling source); 
and the fiducial standard star system. The flux calibrated 
spectrum for a uniform source can be expressed as 

^ C(A) S{X) T{A) 10O-4xr(A) ^ 

in which C(A), is the wavelength calibrated spectrum in 
DN sec“^A“^ per pixel, S'(A) is the sensitivity function 
in ergs DN“^, r(A) is the aperture correction for the slit 
size in question, x is the airmass of the observation, t(A) is 
the extinction correction expressed in mag airmass“^, and 
D is the solid angle of each pixel in steradians per pixel. 
We discuss each component of the calibration separately 
below. 


3.6.1. Point Source Calibration 

We observed Hamuy etal. (1992, hereafter H92) ter¬ 
tiary spectrophotometric standards roughly 15 times each 
night. For each standard observation, we first centered the 
star in a narrow slit (~ 2arcsec), and opened the slit to 
10.8 arcsec after guiding was established. Nine different 
standards were observed throughout the run, with colors 
0.0 < {B — V) < 0.6 mag. This range represents as broad 
a distribution as is available from the H92 standards and 
does bracket the color of the night sky. All standards were 
observed with the slit aligned along the parallactic angle 
(Filippenko 1982). 

One-dimensional spectra were extracted from the two- 
dimensional images and wavelength calibrated in the usual 


way using the APSUM and DISPCOR tasks in IRAF. Ex¬ 
tinction corrections were calculated from standard star 
observations themselves in 50A bins for each night indi¬ 
vidually. Typical residuals in the extinction solution in 
each 50A bin were 0.009 mag {rms) in the data taken on 
1995 November 27 and 29. As no difference was found be¬ 
tween the two nights, the final extinction solution was con¬ 
structed from the data taken on both nights together using 
a sixth-order Legendre polynomial fit to the extinction as 
function of wavelength (see Figure with an uncertainty 
of 0.2% as a function of wavelength. The resulting ex¬ 
tinction curve is in excellent agreement with the r- and 
g-band extinction terms which we obtained from images 
taken simultaneously with the Im Swope telescope at the 
same site. 

After the extinction corrections were applied to the stan¬ 
dard star spectra, the sensitivity curve for each night was 
determined using the task SENSFUNC in the IRAF SPE- 
CRED package. The agreement between the sensitivity 
curves for the two nights is excellent: the variations be¬ 
tween the sensitivity curves found for 1995 November 27 
and 29 are 3% at 4100A, and 2% red-ward of 4500A (see 
Figure Variations in the quantum efficiency of the 
CCD on this level are expected to result from tempera¬ 
ture changes of a few degrees (M. Blouk, personal com¬ 
munication), so we do allow the sensitivity function so¬ 
lutions to be slightly different from night to night. The 
standard deviation in the SENSFUNC solution for both 
nights is 0.011 mag, which translates into a 0.3% error in 
the mean sensitivity as a function of wavelength from 15 
standard star observations. To be conservative, we adopt 
a systematic uncertainty in the point source calibration of 
0.6%. Systematic errors in the tertiary sta ndard star sys¬ 
tem of Hamuy etal. are discussed in §3.7. They are not 
included in the accuracy of the zodiacal light measurement 
discussed here. Those uncertainties are, however, relevant 
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Fig. 6.— The upper plot shows encircled energy curves for stars 
imaged directly in the focal plane of the du Pont 2.5m telescope. The 
half—width of the slit (5.8arcsec) used for standard star observations 
is indicated by the solid vertical line. The percentage flux from a 
point source which passes through this aperture is 97.9 (±0.2)%. 
The lower plot shows the encircled energy curve along the spatial 
extent of a spectroscopic image, showing that 98.4 (±0.1)% of the 
light entering the spectrograph from a point sources is recovered 
when a ±10 arcsec extraction aperture is used. 


to the EBL detection for which we have used this zodia¬ 
cal light measurement. That error is explicitly included in 
final accuracy of the EBL detections discussed in Paper 1. 

3.6.2. Aperture Correction 

The aperture correction compensates for flux which is 
lost from the point source observations in two distinct 
ways. First, light is lost in the focal plane of the tele¬ 
scope if the radius of 100% encircled energy is larger than 
the half-width of the slit. Second, when a one dimensional 
spectrum is extracted from the two dimensional, dispersed 
image, some light will lie outside the extraction aperture 
in the spatial direction. 

To measure flux lost in the aperture plane, we measured 
the PSF from images taken on the last night of the run 
(1995 November 30). We plot the encircled flux with ra¬ 
dius (growth curve) for five stars in Figure^. To be certain 
that focus does not affect the PSF at a 5 arcsec radius, the 
stars used for this plot were taken with widely varying 
focus. As 5 arcsec in the focal plane of the duPont 2.5m 
Telescope corresponds to 0.45 mm (almost 19 pixels) it is 
difficult to have the telescope out of focus enough to affect 
the enclosed flux at a radius of 5 arcsec. This is clear from 
the lack of variation in the shape of the growth curves 
shown in the figure. The fractional flux enclosed by the 
half-width of the slit used for standard star observations 
is 0.979 (±0.002). 

To determine the percentage flux lost once through the 
slit, we have mapped the growth curves along the spatial 
extent of the two dimensional spectra for all 45 standards 
taken during the run. The sky level for this test was taken 
from pixels further than 80 arcsec from the peak of the 


star. We averaged 400 rows near the peak sensitivity of the 
spectrograph to increase the signal-to-noise. To confirm 
that the PSF has negligible wavelength dependence, we 
calculated the aperture correction at both the blue and red 
ends of the spectral range and found no variation in the 
growth curve. The growth curve along the slit for the 18 
highest signal-to-noise spectra are shown in Figure^. The 
differences in the encircled energy at the inner-most radii 
plotted are a result of differences in sub-pixel centering 
for spectra observed with 4x1 binning. This does not 
affect our results for spectra extracted to an aperture of 
±4 (binned) pixels (±10.6 arcsec). As this plot shows, 
the extraction aperture is well within the region of good 
signal-to-noise in the star, and includes 98.4 (±0.1)% of 
the flux from a point source which passed through the 
spectrograph slit. Sky background was measured outside 
an annulus of 25 pixels (60 arcsec) from the peak of the 
star. This aperture is small enough to ensure adequate 
signal-to-noise in the extracted stellar spectrum, and also 
ensures minimal error due to sky subtraction, as the stellar 
signal is much brighter than the sky background in the 
inner 10.6 arcsec. 

The total aperture correction for a uniform surface 
brightness, aperture-filling source is then the fractional 
flux recovered for a point source, or T = 0.963(±0.002). 

3.6.3. Solid Angle of the Program Observations 

The solid angle is a function of both the angular pixel 
scale (spatial direction) and the angular slit-width (disper¬ 
sion direction). The pixel scale was measured empirically 
by taking spectra of two stars with known angular sepa¬ 
ration while they were simultaneously aligned in the slit. 
The measured separation in pixels was then compared to 
the known angular separation of the stars. Four pairs of 
stars were observed in this way with angular separations 
in the range 44-82 arcsec. Each pair of stars was observed 
with the slit at 3 positions differing by less than 1 de¬ 
gree, in an attempt to obtain truly parallel slit alignment. 
The stars used for this purpose are in the field of M67, 
for which the relative astrometry of members is known to 
better than 0.3mas (Girard etal. 1989). The pixel scale 
was found to be 0.5843 ± 0.0035arcsec pixel”^ (Icr error). 

The slit-width of the spectrograph is adjusted by a man¬ 
ual micrometer while the instrument is on the telescope. 
The calibration and repeatability of the micrometer was 
verified using a microscope. The repeatability of the width 
setting was tested by opening the jaws to their maximum 
extent between each of several sets of measurements. At 
the micrometer setting used for our program observations, 
the slit width was measured to be 1.536 ±0.002 arcsec (Icr 
error) at the center of the jaws. The total solid angle of 
each pixel is therefore 0.8975 ± 0.0054arcsec^(lcr error). 

3.7. Accuracy of Tertiary Standards 

Hamuy etal. (1992) quote the internal precision of 
their tertiary spectrophotometric system to be better than 
0.01 mag at all wavelengths, a claim which is well corrob¬ 
orated by the small statistical errors we find in our own 
spectrophotometric flux calibration (see below). An addi¬ 
tional concern for our purposes, however, is the agreement 
between the H92 system and the primary calibration for 
Vega in the spectral range of our observations (approxi¬ 
mately the Cousins B-band, 4200-5100A). 
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ISL) will not suffer the same “extinction”; rather, it will 
appear to be smoothed out over the sky as light from dif¬ 
ferent regions scatters into and out of the line of sight. The 
efficiency of scattering in the atmosphere is conveniently 
described in the familiar way as the extinction along the 
line of sight Tobs(A). 

If the scattering angles were small, and the diffuse, ex¬ 
tended source were uniform over the sky, then the scat¬ 
tering into and out of the line of sight would roughly can¬ 
cel. However, Rayleigh scattering occurs over very broad 
angles and the relative surface brightness of ZL and ISL 
changes strongly over the sky. We must therefore explic¬ 
itly calculate the net effects of scattering in our observing 
situation (determined by the observatory location and po¬ 
sitions of the Sun, Galaxy, and target), as they are not 
intuitive. This is done in detail in the Appendix. We 
summarize the results here. 

We can describe the observed spectrum of the night sky 
in the target field , Ins, as follows: 


Fig. 7. — In the upper plot, the spectrum of the integrated starlight 
(ISL) at |f)| = 60° is compared to the solar spectrum, which has 
been scaled to the same flux and offset to allow visual comparison 
of the spectral features. The lower plot shows the ratio of the two, 
normalized at 4600A. Both plots show visually and quantitatively 
that the absorption lines in the ISL spectrum are weaker than the 
same features in the solar (and therefore zodiacal light) spectrum. 


The tertiary system of H92 is calibrated based on the 
equatorial secondary standards of Taylor (1984), which is 
in turn calibrated to the primary calibration of Vega by 
Hayes & Latham (1975). As described in H92, they recal¬ 
ibrate the Taylor (1984) secondary standards to the now 
widely-accepted Hayes (1985, hereafter H85) calibration of 
Vega. Details of that process are documented in H92. H92 
estimate that their internal consistency in converting Tay¬ 
lor (1984) to the Hayes (1985) Vega system is 0.009 mag in 
the wavelength range of the Cousins B. They also compare 
synthetic photometry of the adjusted Vega spectrum with 
Johnson & Harris (1954) photometry of Vega and find an 
offset of —0.016 (±0.009) mag at B (in the sense of Taylor 
minus Johnson). This offset seems to be intrinsic to the 
Vega calibration, as Hayes also finds that offset between 
his observations of Vega and the original Johnson observa¬ 
tions. We therefore conclude that the statistical accuracy 
of the tertiary system is roughly 1% and the systematic 
uncertainty is roughly 1.5%. 


4. COMPONENTS OF THE NIGHT SKY SPECTRUM 


The spectrum of the night sky from Earth includes inte¬ 
grated starlight (ISL), zodiacal light (ZL), EBL, and dif¬ 
fuse Galactic light (DGL) from above the atmosphere, as 
well as atmospheric emission, or airglow. Light from all 
sources is scattered by molecules and particulates in the at¬ 
mosphere, causing considerable redirection from one line of 
sight to another. Absorption by particulates is a minor ef¬ 
fect compared to the scattering by both components. The 
net effect of atmospheric scattering on a source with small 
angular extent (e.g. stars and galaxies) is atmospheric 
“extinction.” The net scattering out of the line of sight 
can be me asured in the typical way using standard stars, 
as in §3.6. Light from a very extended source (e.g. ZL, 


/Ns(A,t,x) = /(3h.-20d)(A)e-"°-(^)>^ 

±7scat(A, t, x) ± ^air (A,t,x), (3) 

in which /( 3 h,- 20 d) (A) is the flux from the target field (coor¬ 
dinates a = 3.00h, 6 = — 20.18d), Tobs(A) is the extinction 
for a point source, x is the airmass of the target field at the 
time of observations, /scat (A, t, x) is the light scattered into 
the line of sight as a function of time and wavelength, and 
/air(A,t, x) is the effective airglow along the line of sight 
(including any scattering effects which redistribute airglow 
over the sky). The flux from the target field, /( 3 h,- 20 d)(A), 
includes the EBL, DGL and ZL flux within the solid angle 
of the slit. There is no ISL component coming directly 
from the target field because the slit simply contains no 
stars to V = 24 mag. The slit also provides an extremely 
effective pupil stop which prevents contributions from dis¬ 
crete sources off-axis.|^ /scat(A, t,x) can be expressed as 

/scat(A,t, x) = 4^at(A, Xi ZL) ± t, X, ZiL) ± 

(A, t, X, ISL) ± (A, t, X, ISL)(4) 

where the superscript R or M denotes Rayleigh scattering 
(due to molecules) or Mie scattering (due to particulates), 
the parenthetical ZL or ISL denotes the source being scat¬ 
tered, and the parenthetical A, t, and x denote dependence 
on those variables. We do not include the DGL in Equa¬ 
tion ^ explicitly because the total DGL is at least a factor 
of 50 times fainter than the direct ISL and is therefore a 
trivial component (< 0.2% at 2<UT<6.5) of the scattered 
light. The EBL is not included as it does not have strong 
spectral features (see §||). 

At any altitude, Rayliegh scattering is the dominant ef¬ 
fect. Because the particulate density is concentrated at low 
altitudes, this is especially true at high altitude observato¬ 
ries. The total extinction, as measured for a point source, 

^No stars with U < 12 mag are within 12 arcmin of the slit, 
and no stars with V < 7 are within 1.5 degrees. We have carefully 
characterized the scattered light properties of the duPont telescope 
by positioning a U = 4 mag star around the field from on-axis to 
20 arcmin off-axis in 4 directions at 1 arcmin intervals. The stray 
light entering the slit from discrete off-axis sources is more than 10“® 
fainter than surface brightness of the ZL in the field. See Paper I for 
further discussion. 
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is equal to the sum of the molecular and particulate extinc¬ 
tion, Tobs(A) = 'rM(A)-l-Tp. (A). Rayleigh extinction, tr, can 
be calculated from the well-known density distribution of 
the atmosphere for any observatory. Mie extinction, which 
varies with time and geography, can be inferred from the 
difference between the observed and Rayleigh extinction. 
At LCO, the extinction due to Mie scattering is 20-40% 
of the Rayleigh extinction. 

We have calculated the scattered light from all terms 
in Equation ^ in the Appendix at 30 minute intervals 
throughout the nights of our observations. To briefly sum¬ 
marize the results of our calculations, the ZL total flux 
scattered into the line of sight {gained) at any time dur¬ 
ing our observations is less than the total flux scattered 
out of the line of sight {lost). Thus, the net result of at¬ 
mospheric scattering for the case of ZL in our situation is 
still a net extinction of order 2-8%, which we can conve¬ 
niently describe by an effective extin ction, r^ff , whi ch we 
use in place of Tobs for ZL (see Figures AlO and All in the 
Appendix). We can check the scattering predictions of our 
calculations in our ZL analysis itself by looking for changes 
in the ZL solution with time. We estimate that our cal¬ 
culation of the scattered ZL has an average uncertainty of 
8%, which translates into a systematic uncertainty in our 
ZL measurement of 1.2%. 

In the case of the ISL, the total flux gained due to scat¬ 
tering into the line of sight is 12-24% of the total ZL flux 
from low to high airmass. However, the crucial issue is 
not the total mean flux, but rather the strength of the 
spectra l feat ures which are in common with the Sun (see 
Figure A15 in the Appendix). The net influence of the 
scattered ISL on observations is to increase the strength 
of the Fraunhofer lines over the night by 0.6 — 4% red- 
ward of 4I00A, and 5-35% blueward. Because the effect 
is a strong function of wavelength, it is straightforward to 
identify inconsistencies between the predicted ISL flux and 
our observations by looking for changes in the ZL solution 
with wavelength. We estimate that our calculation of the 
scattered ISL has an uncertainty of 13%, which translates 
into a systematic uncertainty in our ZL measurement of 
0.5% over the majority of our wavelength range. 


5. ANALYSIS 

Based on the discussion in the previous section, we can 
characterize the observed spectrum of the night sky, Jns, 
at any time as a combination of airglow, zodiacal light 
(ZL), diffuse Galactic light (DGL), and scattered inter¬ 
stellar light (ISL) as follows: 


^Ns(A,t, x) — ^air (A,t,x) -k I scat (A,t,x,ISL) 
+/zL(A)e-"^«(^’*)^ 

+ [.fEBL(A) -k/DGL(A)]e“'^°*’'‘^'^^^, (5) 

in which t^(A, t) is now the effective extinction for ZL dis¬ 
cussed in and /scat (A, t, x, ISL) is the total scattered ISL 
flux due to Rayleigh and Mie scattering in the atmosphere. 

Sharp spectral features are not expected in the EBL be¬ 
cause redshifting will blur any distinct spectral features. 
The EBL will therefore not affect our measurements of 
the ZL from the strength of the observed Fraunhofer lines. 
The diffuse Galactic light, which results from scattering of 
the ambient interstellar radiation field by interstellar dust, 


is very weak in our field (0.8% of the ZL flux, see Paper 
I). In addition, as we have already discussed regarding the 
scattered ISL, the strength of the spectral features we use 
in our analysis is roughly 1.5 to 3.8 times weaker in the 
DGL than in the ZL spectrum. The DGL from the target 
field thus contributes at most 0.2-0.5% percent to the final 
result (see Figure]^). We subtract this contribution from 
our ZL measurement after the fact at a level of 0.3%. Emis¬ 
sion lines due to ionized gas in the DGL do not contribute 
in the spectral range of these observations (see Paper I, 

Martin etal. 1991, Dube etal. 1979, and Reynolds 1990). 

The observed spectrum can therefore be expressed as 
the sum of four components: (a) an unstable emission 
line spectrum, due to airglow; (b) a stable and feature¬ 
less component, due to EBL; and (c) a stable, absorption 
line component, due to ZL; and (d) a time variable ab¬ 
sorption component due to scattered ISL. The component 
(c) can be ignored, and (d) has been calculated. The por¬ 
tion of the night sky spectrum which has variable spectral 
features can therefore be expressed as 

/obs(A, t, x) = /air(A, t, x)+/scat (A, t, X, ISL)+c(A)/q (A)e-"»«fo’‘)^ 

( 6 ) 

in which (A) is the solar spectrum and c(A) is a scaling 
factor which relates the mean surface brightness of the ZL 
to the mean flux of the Sun. 

Identifying the appropriate scaling spectrum, c(A), is 
complicated by the fact that the ZL — a pure absorption 
line spectrum — is greatly obscured by the emission line 
spectrum of the airglow in the night sky. Airglow features 
do overlap with some of the solar Fraunhofer lines, as can 
be seen in the comparison of a night sW spectrum and a 
scaled solar spectrum shown in Figure g. The strength of 
particular airglow lines can vary by several percent during 
a single night, as can be seen in the comparison of two 
night sky spectra shown in Figure |^. The airglow spec¬ 
trum is composed of an effective continuum due to O-kNO 
(NO 2 ) recombination, broad rotation-vibration transition 
bands, scattered light, and blended lines, making a con¬ 
tinuum level impossible to identify. Not only do rapid 
temporal variations occur with (de)ionization, but airglow 
is also a complex function of airmass, observatory location, 
and local atmospheric conditions such as volcanic activity 
(van Rhijn 1924, 1925; Roach & Meinel 1955). In short, 
a stable, fiducial airglow spectrum with meaningful abso¬ 
lute or relative flux cannot be defined. As a result, it is 
not possible to measure the equivalent width of individual 
ZL Fraunhofer lines because the ZL continuum level is well 
hidden. Instead, we have developed a conceptually simple 
approach to the problem of determining the scaling fac¬ 
tors c(A) which does not involve measuring the equivalent 
widths explicitly. 

We begin with the assumption that the intrinsic airglow 
spectrum, time- and airmass-dependent though it may 
be, does not have spectral features in common with the ZL 
spectrum. This is borne out by the fact that we get consis¬ 
tent solutions usiM eight distinct spectral regions spread 
over llOOA (see (Q. When we subtract /scat(A,t,X 5 lSL) 
and scaled solar spectrum with the correct value for c(A), 
what remains is a pure airglow spectrum, free of solar fea¬ 
tures: 


/air(A,t,x) — Aibs(A,t,x) Acat (A, t, Xj I^L) 
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Wavelength (A) 

Fig. 8.— The observed spectrum of the night sky compared to a solar spectrum at arbitrary absolute flux. The solar Fraunhofer absorption 
lines which are preserved in the ZL are clearly visible in the spectrum of the night sky. The strongest of these lines are marked; several blended 
solar features are also seen in the night sky spectrum. Square brackets indicate the wavelength regions used in the final analysis. 


-c(A)/Q(A)e-"-(^’‘)>^. (7) 

We use a linear correlation function to determine when 
the difference (residual airglow) spectrum is uncorrelated 
with the solar spectrum and is consequently free of so¬ 
lar features. When the correlation between the difference 
spectrum and solar spectrum is minimized, the correct ZL 
surface brightness has been subtracted from the observed 
night sky. 

The only available, high-resolution spectrum of the Sun 
is the National Solar Observatory Solar Flux Atlas of the 
integrated solar disk at O.OlA resolution. The statistical 
error in the flux calibration of this spectrum is 0.25% as 
estimated by agreement in overlapping sections of the nor¬ 
malized spectrum. The effects of atmospheric absorption 
by H 2 O or O 2 are negligible below 6500A, as described 
in the published Atlas (Kurucz etal. 1984). In the opti¬ 
cal, the normalized spectrum can be converted to absolute 
solar irradiance using the Neckel & Labs (1984) (NL84) 
absolute calibration. As the NL84 is the standard with 
respect to which the ZL color is defined, the absolute ac¬ 
curacy of the fiducial solar spectrum is not a source of error 
in this work. The Solar Flux Atlas, calibrated to NL84, 
was obtained in digitized form from R.L. Kurucz. It was 
convolved with a variable width Gaussian to match the 


resolution of the observed spectra as a function of wave¬ 
length. The wavelength-dependent resolution of each pro¬ 
gram spectrum was determined from the arc lamp spectra 
which were used for wavelength calibration. 

The execution of this method is complicated by the fact 
that the relative color of the airglow and solar spectra 
will dominate the strength of the diagnostic correlation 
if the continuum shapes of both spectra are not properly 
removed. In order for the strength of the linear correlation 
of /air (A) and Iq (A) to reflect the strength of coincident 
spectral features, both spectra must have stationary mean 
values as a function of wavelength, as can be seen clearly 
in the generic expression for a linear correlation: 


R{x,y) 


-x){yn-y) 

\/Y.ni^n - ^) VE„(j/n “ V)^ 


( 8 ) 


In this case, x and y are /air(A) and Iq (A) respectively, 
the subscript n runs over wavelength. It is clear from this 
expression that the mean flux drops out of the correlation, 
while differences from the mean are crucial. 

Of the 47 strongest solar features, we find that 39 give 
ZL solutions which vary with time by more than 18% over 
the night. In all cases, this variation is correlated with the 
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Wavelength (A) 

Fig. 9.— Two spectra of the night sky taken on the same night, several hours from twilight. Rapid fluctuations are evident in the strength 
of many of the airglow features. 


strength of adjacent airglow lines. The remaining eight so- 
lar features vary by less than 10% with time. The results 
discussed below are based on these eight spectral regions, 
indicated in Figures || and In these regions, the contin¬ 
uum can be well approximated by a simple second order 
polynomial fit and easily subtracted, however our results 
are quite insensitive to the method of continuum fitting. 
Boxcar smoothing with scales between 75A and 199A (at 
least twice the width of the widest spectral region used in 
the analysis), second or third order polynomial fitting, and 
Savitsky-Golay smoothing (Press et al. 1992) all produce 
identical results. 

6. RESULTS 

In Figure |^, we show an example of the correlation 
strength, R, and correlation probability, P, between a ZL 
spectrum and the spectra which result when we assume a 
range of values (8xl0“® to 1.4x10“^ ergs s“^ cm“^ sr“^ 
A“^) for the mean ZL flux contributing to one region of 
one observed night sky spectrum. Where R and P go to 
zero, the correct ZL surface brightness has been removed 
from the observed night sky spectrum. 

The absolute flux of the ZL as measured from each of 
the 16 spectra taken on 1995 November 27 and 29 are 
shown in Figures In Figure |l|, we show the av- 



8x10-® 10-’ 1.2x10-’ 1.4x10-' 

IjL (ergs/sec/cmVst/A) 

Fig. 10.— The value of the correlation parameters used to de¬ 
fine the strength of the correlation between the residual “airglow” 
spectrum (the zodiacal light-subtracted night sky spectrum) and the 
ZL spectrum for different assumed contributions of ZL. The points 
connected by the dotted line indicate the probability that the two 
spectra are from the same parent set as a function of the ZL flux 
assumed. The points connected by the solid line indicate the corre¬ 
lation strength, zero being no correlation. These parameters describe 
a simple linear correlation as defined in equation N (see Press et al. 
1992). The ZL surface brightness identified by this method for the 
observation shown here is 1.08x10“^ ergs s“^ cm“^ sr~^ A“^. 


erage c(A) as determined for all 16 spectra over the run 
for each of the eight spectral regions. The results are nor¬ 
malized to 1.0 at 4650A. This plot illustrates two points: 







































12 


Detections of the Optical EBL: Measurements of Zodiacal Light 



Fig. 11.— The joints show c(A) for the eight spectral regions 
indicated in Figure g as measured from 16 spectra taken on 1995 
November 27 and 29, normalized to 1.0 at 4600-4700A. The error 
bar on each point indicates the standard deviation in the 16 mea¬ 
surements. The hatched region shows the color of the ZL relative 
to the solar spectrum determined by a linear least squares fit to 
the points, excluding the point at ~ 4700A. The best fit color is 
C(3900,5100) = 1.05 zb 0.01, where the error corresponds to the 
one-sigma error in the fitted slope. The mean flux of the ZL at 
4600-4700A from this fit is 109.1(di0.5) x 10“® ergs s“^ cm“^ sr“^ 

A-L 



Airmoss 


Fig. 12.— The mean value of the ZL at 4600-4700A as measured in 
each of the 16 spectra taken on 1995 November 27 and 29, adopting 
C(3900, 5100) = 1.05. Each point represent the average of the ZL 
flux measured in all 8 spectral features in a single observation. The 
horizontal line shows the mean ZL flux, which is 109.7(±0.7) X 10~® 
ergs s“^ cm“^ sr“^ A“^. Dashed lines show the one-sigma statistical 
error in the mean (0.6%). 


(1) from any single spectral feature, we find a solution for 
the ZL flux with a standard deviation of only 1-6% over 
16 individual observations taken on two nights and inde¬ 
pendently reduced; and (2), the solution from independent 
spectral features are in excellent agreement, and indicate 
that the ZL is roughly 5 ± 1% redder than the solar spec¬ 
trum, or (7(3900, 5100) = 1.05 ± 1 (excluding the point at 
~ 4700A). This color is in excellent agreement with pre¬ 
vious estimates from the ground and space (see Leinert 
1998). It is also in excellent agreement with our own mea¬ 
surement of the ZL color from our simultaneous HST/FOS 
observations, which give (7(4000,7000) = 1.044 ± 1 (see 
Paper I). 

From the solution for the ZL as a function of wave¬ 


length we can explore the impact of the scattered ISL flux 
on our measurements. If we increase (decrease) the to¬ 
tal flux of the predicted scattered ISL, /scat(A, t, x, ISL), 
the solution for the ZL decreases (increases) linearly in 
response; a change of 50% in the ISL flux corresponds 
to a change of 6% in the ZL solution at ^ 3950A, but 
< 2% at the other wavelengths. Thus, increasing or de¬ 
creasing /scat (A, t, X, ISL) by 50% makes the ZL solution at 
3950A inconsistent at the two-sigma level with solutions 
over the rest of the spectrum for a ZL color of 5 ± 1%. 
Also, increasing or decreasing the ISL flux consistently in¬ 
creases the scatter in the ZL solution at all wavelengths; 
at 3950A, the scatter increases by 40% in response to a 
change of 50% in the ISL flux at all airmasses. Although 
this does not allow us to place a stronger constraint on the 
error in /scat(A, t, Xj ISL) than those discussed in the Ap¬ 
pendix, it does provide independent verification that the 
predicted /scat (A, t, Xi ISL) values are in the right range. It 
also emphasizes that an error of 10% in /scat(A,t,X 5 lSL) 
changes the mean ZL solution by only 0.4% at 4200- 
5100A. As discussed in the Appendix, the uncertainty in 
/scat (A, t, X) ISL) contributes an uncertainty to the ZL mea¬ 
surement in §^f < 0.5% (see Table ||). 

In Figures O and 13 we show the mean ZL solution 
at 4600-4700A(for (7^00, 5100) = 1.05) as a function of 
airmass from each of the 16 exposures. Figure ^ shows the 
solution obtained using all eight spectral features solutions 
plotted in Figure |^, while Figure ^ shows the solution 
based on the four spectral features with the smallest stan¬ 
dard deviations in Figure |^. The horizontal line in each 
plot shows the mean, the dashed lines show the one-sigma 
error in the mean. The difference between the results in 
the two plots is less than 0.2%. No obvious trends appear 
in either plot as a function of airmass. In fact, the mean 
value for points above and below 1.2 airmasses in Figure |l3| 
agrees to better than 0.3% (Lalf of the error in the meanp 
The error bars in Figure O vary from 1-8%, indicative 
of the small number of measurements (four spectral fea¬ 
tures) being averaged together to produce the result at 
each airmass. The standard deviations indicated by the 
error bars in Figure ^ show less variation from point to 
point (1.5-5%), as eight measurements contribute to each 
point. 

As discussed in the Appendix, we estimate that the un¬ 
certainty in the calculated scattered ZL flux contributes an 
uncertainty to the ZL solution of 1.2%. The stability of 
our ZL solution with airmass indicates that our calculated 
net extinction, which incorporates the ZL scattering mod¬ 
els from the Appendix, has the correct behavior over the 
night. However there is no way to independently infer the 
accuracy of the mean net extinction from the rms scatter 
in the ZL solution because the spectral shape of Teff(A,t) 
does not change significantly over the night (see Figure 
All); an error in the mean level of Teg will not affect the 
rms scatter in the solution between exposures. The errors 
in our result are summarized in Table |^. We discuss the 
systematic uncertainties further in the next section. 


7. DISCUSSION 

The rms scatter in our ZL solution is less than 1%. 
This demonstrates that statistical errors are quite small, 
be they a result of instrumental effects or our analysis 
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Wavelength (A) 

Fig. 14.— Airglow emission lines in the observed spectrum of the night sky after zodiacal light is subtracted. Broad emission features 
(marked above) are molecular rotation—vibration bands of O 2 , N 2 , H 2 , OH, and N02- Collisional de-excitation of these molecules contributes 
a continuum as well. Emission lines (marked below) are also seen in this wavelength range from atomic transitions (photoionized O and Hg). 
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Fig. 13.— The same as Figure hd, but here only the four features 
with the smallest error bars havebeen used to produce the mean 
value of the ZL at 4600—4700A in each spectrum The horizontal line 
shows the mean ZL flux, which is 109.4(±0.6) X 10“® ergs s“^ cm~^ 
sr“^ A“^. Dashed lines show the one-sigma statistical error in the 
mean (0.6%). 


method. Systematic uncertainties due to instrumental ef¬ 


fects are also quite small and are straightforward to quan¬ 
tify (see Table H). This is demonstrated by the fact that we 
obtain consistent results to within 0.3% from data taken 
on two different nights, independently reduced and cal¬ 
ibrated. The uncertainties in the atmospheric scatter¬ 
ing model described in the Appendix have been consid¬ 
ered very carefully and we believe that the adopted un¬ 
certainties are conservative. However, the measurement 
presented here is obviously complex and might be effected 
by systematic errors which are more difficult to anticipate 
or quantify. 

One such systematic effect might include moonlight or 
sunlight scattered in the atmosphere. Based on the scat¬ 
tering analysis in the Appendix, it is clear that scattering 
into the line of sight from any source, even the sun or 
moon, is negligible when that source is more than 14 de¬ 
grees below the horizon. All of the observations in this 
work took place when the sun was more than 18 degrees 
below the horizon. Although the observations took place 
several days after new moon, the moon was below the hori¬ 
zon during all of our program observations, and below 14 
degrees for all but 1 exposure. In addition, the net effect of 
such solar-type scattering contributions, if present, would 
be to increase our estimate of the zodiacal light, and con- 
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sequentally to artificially decrease the value of the inferred 
EBL in Paper 1. We do not believe that such scattering is 
likely to have influenced our results. 

As discussed in §§|| and |[ the contamination of solar 
features by airglow, while introducing a systematic error, 
would not introduce a stable systematic error: the flux 
of airglow features changes constantly through the night. 
The stability of our ZL solution in the eight spectral re¬ 
gions we have used demonstrates empirically that airglow 
is unlikely to have had a significant effect on our results. 
However, the possibility can’t be ruled out and may intro¬ 
duce a systematic error which we cannot quantify and is 
not included in our estimate of the formal errors. 

It is also possible that some Doppler shifting occurs in 
the ZL spectral features relative to the solar spectrum due 
to motion of the dust in the zodiacal plane. Eor that 
reason, we allowed for a shift in the central wavelength 
when calculating the correlation but found no measurable 
offset. We note, also, that the results of this method 
would not be affected by the slight Doppler broadening 
which might affect the spectral features of the ZL, because 
Doppler broadening will not alter the total flux across a 
feature. The correlation is unaffected by the saw-tooth ef¬ 
fect of subtracting features with mismatched widths at the 
level of the O.SA Doppler broadening which is expected at 
the orientation of these observations (East & Reay 1984). 
Note also that while the resolution of the input spectra 
used for calculating the scattered ISL flux is lower than the 
resolution of our program observations (4A versus 2.6A), 
this will not affect our analysis as long as regions with 
width >>4A are used in the analysis. The smallest of our 
spectral regions is ISA. 

Figure H shows the airglow spectrum (the night sky 
spectrum after zodiacal light is subtracted) we obtain 
by this method. Emission lines from molecular rotation- 
vibration bands (O 2 , N 2 , H 2 , OH, and NO 2 ) are labeled, as 
are some atomic transmission lines (O and Hg). Identifi¬ 
cation of emission features in this range of the spectrum is 
not complete (see Schmidtke etal. 1985, Slanger & Huestis 
1981, Jones etal. 1985 and references therein). 

Finally, we note that our measurement of both the mean 
flux and ZL color are in very good agreement with typical 
values for the similar viewing geometries quoted in the 
literature (see the results of Levasseur-Regourd & Dumont 
1980, pictured in Figure A4, and Leinert 1998). 


8. SUMMARY 


We have measured the mean surface brightness of the 
zodiacal light along a single line of sight towards an ex- 
tragalactic target using ground-based spectrophotometry 
with a 300 arcsec^ field of view. The observations were 
made on on 27 and 29 November 1995, simultaneous with 
HST observations of the same field. The goal of this co¬ 
ordinated program is a measurement of the optical extra- 
galactic background light. Because the zodiacal light at 
optical wavelengths results from sunlight scattered by in¬ 
terplanetary dust concentrated in the ecliptic plane, the 
flux towards an extragalactic field varies seasonally. Vari¬ 
ations in the interplanetary dust cloud with time and the 
solar cycle may also affect the flux of zodiacal light. The 
measurement of zodiacal light present here is therefore 
uniquely relevant to the date and target field. 

Our results incorporate explicit calculation of net effect 


of atmospheric scattering on terrestrial measurements of 
the zodiacal light, and show that these effects are small 
(< 10%) for zodiacal light measurements far from the Sun. 
We find the mean flux to be 109.4x10“® ergs s“^ cm“® 
sr“^ A-i at 4650A (see Figure |l^), and the color to be 
5(±1)% redder than the solar spectrum per lOOOA. The 
statistical uncertainty in the mean flux is 0.6% (lu), and 
the systematic uncertainty is 1.1% (Icr). We discuss addi¬ 
tional systematic effects which might influence this mea¬ 
surement beyond those which are quantified here. Our re¬ 
sults are in good agreement with previous measurements 
of the ZL at similar orientations with respect to the eclip¬ 
tic plane and scattering geometry (see Leinert etal. 1998 
for a recent review). This is the only optical measurement 
to date which isolates the ZL from other uniform back¬ 
grounds, including diffuse Galactic light and extragalactic 
background light. The color of the ZL as a function of 
the line of sight through the interplanetary dust cloud is 
further addressed in Paper 1. 
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APPENDIX 

ATMOSPHERIC SCATTERING 

In this Appendix, we present all scattering calculations 
and results used in Section Scattering in the atmosphere 
is a combination of Rayleigh scattering by molecules and 
Mie scattering by particulates. The scattering due to these 
two components can be dealt with individually. We begin 
by describing a general model for the scattering in a spher¬ 
ical atmosphere and then discuss the specific parameters 
needed to calculate Rayleigh and Mie scattering affecting 
observations from Las Campanas Observatory. Finally, we 
present the predicted contribution of scattered light to the 
program observations analyzed in this paper. We address 
zodiacal light and the integrated starlight as scattering 
sources separately. 


Generic Calculations 

The first calculations of radiative transfer in a Rayleigh 
scattering atmosphere were published by Chandresekar in 
1950. Since then, a number of authors have published ra- 
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Observer’s 


Top of the line of sight 

Atmospliere Zenith Zenith (A,Q^={Aq,Zq) 



Fig. A1.— Geometry for atmospheric scattering along the line 
of sight OXT^ with the observer at O and the line of sight exiting 
the atmosphere at T. In the case illustrated, scattering occurs at X 
from light entering the atmosphere at N from N'. 



Fig. A3.— Correction factor for multiple scattering, Fms, as a 
function of the Rayleigh extinction. 



Wavelength (A) 


Fig. A2.— We compare the observed extinction, t(A), derived 
from observations of standard stars on 1995 November 27 and 29, 
with the predicted extinction due to Rayleigh scattering, which is 
well known (see text). The extinction due to Mie scattering is the 
difference between the two. 


diative transfer calculations addressing Rayleigh and Mie 
scattering in a curved atmosphere (e.g. Sekera 1952, Sek- 
era & Ashburn 1953, and Ashburn 1954), and the effects of 
multiple-scattering (Dave 1964, de Bary & Bullrich 1964, 
and de Bary 1964). Careful measurements of zodiacal light 
over the sky and intensity distributions of the daytime sky 
have empirically demonstrated the accuracy of those cal¬ 
culations (e.g. Elterman 1966; Green, Deepak, & Lipofsky 
1971; Weinberg 1964, Dumont 1965). 

To calculate the atmospheric scattering affecting the ob¬ 
servations described in this Paper, we begin by adopt¬ 
ing the scattering geometry and coordinate system defi¬ 
nitions used by Wolstencroft & van Breda (1967, hereafter 


WvB67), illustrated in Figure [aI]: (A, ^)o and {A',4’)x 
are azimuth/zenith-distance coordinate systems centered 
on the observer at O and a generic point, V, along the 
line of sight, respectively. The problem is then to calcu¬ 
late the brightness observed at O, along the line of sight 
("4, C)o = (Ao, zo)o- 

Following WvB67, scattering occurs at the point X for 
radiation which entered the atmosphere at the point N 
from the direction N', given by (A, (C)o or {A',(j))x- The 
light arriving at X from N' can be expressed as 

Ix = L{A\(t>) sincf) (Al) 


in which the above-the-atmosphere source has flux A (A', cj)), 
light is attenuated by g^g j,. ^;|.g^ygig along 

NX, and s is the distance along that path. Attenuation 
is a function of C'ext(A), the extinction cross section of the 
scattering particles in cm^, and of hi{cj),a), the effective 
column density of particles along the line of sight. The ef¬ 
fective column density is defined by the local zenith angle, 
(j), and the distance, a, which defines t he point X relative 
to the center of the Earth (see Figure [aI] ): 

^'Smax ) 

hi{(j),a) = / n{a')ds', (A2) 

Jo 

in which Smax(<(>, o’) is the distance from X to the top of 
the atmosphere at N, and n(a') is the atmospheric num¬ 
ber density of molecules in cm“^ as a function of distance 
from the center of the Earth, a', and as a function of the 
distance s' along the line XN. 

The light scattered towards the observer from X is then 

Gcat(A)n(a)P(d) Jx#dA'ds, (A3) 


in which P{0) is the scattering phase function and Ix^ 
the flux arriving at point X is given in Equation Al. Fi¬ 
nally, the scattered light is further attenuated by the factor 












16 


Detections of the Optical EBL: Measurements of Zodiacal Light 



500 ‘000 1500 


Fig. A4.— Zodiacal light over the sky is plotted in sun-centered 
ecliptic coordinates, so that the Sun is at (A — \q = 0°,/3 = 0°), 
the anti-solar direction is (A — Aq = 180°,/3 = 0°), and the ecliptic 
poles are at /3 = ±90°. All four quadrants are assumed to be the 
same, which is correct to within a few percent during the time of 
year when our observations occur. 


g C'ext(A)^ 2 ( 20 ,O-)^ which 


h2{zo,cr) = hi{zo,cro) - hi{zo - (A4) 


The total flux scattered into the line of sight (zq, Aq) from 
sources distributed over entire visible hemisphere of the 
sky is then 


'fscat('2'0 7 Aq) 


3Cscat(A) 

IGtt 


{zo,(To) ^7r/2 + /(cT) ^2-77 


/s=0 




fA'^0 


n{a) P{0) L{A\ (j)) sincj) (^5) 

g-C.cat(A) [h,{<l>,a) + h2(zo,^)] 


The visible sky at the point X dips below the observer’s 
horizon at large values of s. Hence, the limit of the in¬ 
tegral over (j), is greater than (j)/2 by the value /(cr) = 
cos“^(i?/cr) < 14°, where R is the radius of the Earth 
(6371 km). The equations needed to change variables be¬ 
tween {A',(j)) and (zq, Aq) are given WvB67. 

The phase function for Rayleigh scattering is P{9) = 
1 -|-cos^(0). The atmospheric density is given by the stan¬ 
dard barometric formula n{a) = , where H is the 

altitude above sea level {H = a — R), the scale height is 
Ho = 7.99 km, the density at sea level is no = 2.67 x 10^® 
cm“^ , and the effective scattering cross-section for air 
is Cscat = 7.78 X 10“^'^(A/4600A)“"^cm^ (see Schubert & 
Walterscheid 1999 and van de Hulst 1952). Eor molecules 
in the atmosphere, extinction is entirely due to scatter¬ 
ing, so that Cext = Gcat- Atmospheric extinction due to 
Rayleigh scattering is then r^(A) = Cext(A) n{a)da. 
Eor the duPont telescope at Las Campanas, which is at 
an altitude of 2.28 km, the expected Rayleigh extinction 
is r^(4600A) = 0.12. 

The phase function for Mie scattering by particulates 
in the atmosphere depends on the distribution of parti¬ 
cle sizes, and must be empirically determined. We adopt 
the phase function measured by Green, Deepak & Lipof- 
sky (1971) from their complete analysis of the Mie (par¬ 
ticulate) scattering and Rayleigh (molecular) scattering 
components of the atmosphere based on the scattering of 


sunlight. Their results are in good agreement with theo¬ 
retical scattering models and other estimates of the size- 
distribution of particles in the atmosphere and have a neg¬ 
ligible dependence on wavelength for our purposes (see El- 
terman 1966, and Deepak & Green 1970). The scattering 
and extinction coefficients for particulate scattering are a 
function of the size distribution of particles and vary with 
time and geography. The extinction due to Mie scatter¬ 
ing can, however, be inferred from the observed extinction 
for a point source and the calculated Rayleigh extinction 
coefficient: = Tobs — ~ 0.05 at 4500A (see Figure 

|A2D . This vafue is in good agreement with estimates for 
Tenerife by Dumont (1965) and for Hafeakala by Weinberg 
(1964). This is not surprising as our observed Tobs(A) is 
consistent with the GTIO curves (Bafdwin & Stone 1984, 
Stone & Baldwin 1983), and t^(A) is simply a function of 
the atmospheric density. 

Unlike the case for molecules, the attenuation caused 
by particulates is not entirely due to scattering. Staude 
(1975) adopts values of Cscat = 4.47 x 10“^^ cm^ and 
Cext = 7.53 X 10“^^ cm^ for dry air at 4200A. With a sea 
level density of no = 1.11 x 10^ cm^, and a distribution 
scale height of only ho = 1.2 km, these parameters give 
~ 0.01 at 2 km. We have scaled Cscat and Cext to give 
values consistent with our observed value of r^. Scaling 
Ho or assuming a different value of H would have the same 
effect on J^^A-^^ZL). 

The scattering modef discussed above describes a sin¬ 
gle scattering event. However, multiple scattering events 
become significant for scattering angles 9 > 30° (de Bary 
1964, de Bary & Bullrich 1964). Consequently, we ap¬ 
ply a multiple scattering correction for Rayleigh scattering 
which is adopted from Dave (1964) and plotted in Figure 
A3. The correction factor plotted in Figure A3 is sim¬ 
ply the factor by which the intensity of scattered light 
increases over the single-scattering case. Multiple scatter¬ 
ing does not occur due to particulates (Mie case) because 
of the small scattering angles which dominate that process 
and very small values of (X). 

To confirm the accuracy of our calculations, we checked 
our scattering model against published results of Staude 
(1975), WvB67, and Ashburn (1954) for a uniform, sky¬ 
filling source of unit flux. We find that our results are 
consistent to within 4% for zenith angles z < 80° before 
the multiple scattering correction is applied. (WvB67 pre¬ 
dates evidence for the effects of multiple scattering, and 
Staude adopts the same corrections used here.) The un¬ 
certainty in the multiple scattering correction is roughly 
4-7%, increasing with larger values of r. 

Using the expressions above, we calculate the scattered 
light flux, /scat(A), resulting from Mie scattering by par¬ 
ticulates and Rayleigh scattering by molecules throughout 
the nights of our observations. The results depend explic¬ 
itly on the absolute position of the Sun and the Galactic 
center relative to the observatory and relative to the tar¬ 
get field. In the following sections, we consider the cases 
of zodiacal light (ZL) and integrated starlight (ISL) as the 
extra-terrestrial source of flux separately. 


Zodiacal Light 

To calculate the scattered zodiacal light along the line 
of sight of our observations, we adopt values for the zo¬ 
diacal flux given in Leinert etal. (1998), which are taken 
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UT= 7.5 


Fig. A5.— Integration pattern for calculating scattered zodiacal light. Each plot is an Aitoff projection of the sky in Ecliptic coordinates. 
The position of the Sun and the ecliptic plane are indicated by the Sun symbol with a bar through it. Obviously, the ecliptic plane is along 
the 0° latitude line. Our target field is indicated by the square. The ecliptic coordinates of the zenith can be seen as the “bullseye” center 
of the integration pattern. The motion of the target relative to the local zenith is clear as the target moves from east to west of the zenith 
point in the integration pattern. The coordinates of the local horizon are traced by the edge of the integration pattern. The Sun is obviously 
located just to the west of the horizon at the start of the night (UT=2.0 hr) and just to the east of the horizon at the end of the observing 
night (UT=7.5 hr). 




Fig. A6. — The Rayleigh scattered zodiacal light contribution to 
our observed night sky flux on the nights of 1995 November 27-29 
as a function of airmass at lOOA intervals from 3900A to 5100A. 
Flux decreases monatonically with wavelength at any airmass, fol¬ 
lowing the behavior of r^(A). 7^^^(A,ZL) is given as a fraction of 
the above-the-atmosphere zodiacal light flux in our target field; this 
removes the spectrum of the zodiacal light itself and emphasizes the 
wavelength dependence of the scattering. The each line shows the 
change in ZL) as the target field goes from slightly east of 

zenith, through zenith, to far west of zenith during the night. Dots 
mark half-hour intervals between UT=2.0 hr and UT=7.5 hr. 


Fig. A7. — Same as Figure |Aq , but here we plot the contribution of 
Mie scattered zodiacal light along the line of sight. Again, the scat¬ 
tered light flux is plotted as a fraction of the above-the-atmosphere 
zodiacal light flux in our target field. 


from Levasseur-Regourd & Dumont (1980) with values at 
elongations e < 30° added from space-based observations. 
These ZL values are above-the-atmosphere fluxes and are 
in excellent agreement with later space-based results (see 
Leinert etal. 1981). To obtain a s mooth flux distribution 
of ZL on the sky (see Figure A4), we use the spherical 
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Fig. A8.— Same as Figure but here each line indicates 
^iSat ^ function of wavelength at discrete times. Scattering 

is maximized at high airmass (far west of zenith), and minimized 
slightly west of zenith, when the airmass is still relatively low and 
the Sun has had time to set far enough that highest flux regions of 
zodiacal light are no longer in the visible hemisphere of the sky. 



Wavelength (A) 

Fig. A9.— Same as Figure |A^ , but showing the Mie scattered 
zodiacal light. 


interpolation me thod developed by Renka (1997). 

In Figure we show the integration pattern in cj) 
(zenith an gle) and A (azimuth) used to calculate /scat('^) ZL) 
(Equation |A6[ ) at the indicated times during the nights of 
our observations. Actual calculations were done with twice 
the number of integration points shown in the figures. The 
visible part of the sky (shown by the integration pattern) 
is at least 30 degrees from the Sun during our observations. 

As a technical detail, we have made the simplifying as¬ 
sumption that the spectral shape of the ZL over the visible 
hemisphere is uniform. That is, only the mean flux of the 
ZL changes. Although there are variations in the color (de¬ 


fined in Equation]^) of the ZL over the range 3900-5100A 
from e = 30° to e = 180°, the total change is empiri¬ 
cally less than ~ 8% and our target is in the center of the 
expected color range (e.g., Frey etal. 1974, Leinert etal. 
1981). We have run trail scattering models in which we 
change the flux with e over the sky by ±4%, and we find 
that the effect on the predicted scattered flux is changed 
by 0.4% at airmasses higher than 1.6, and 0.2% at the 
lowest airmass. In other words, by ignoring the color vari¬ 
ation in ZL over the sky, the scattered light model will be 
wrong by 0.2% at 3900A relative to the value at 5100A, or 
±0.1% over the range 3900-5100A for our observing situ¬ 
ation (positions of the Sun relative to the target and the 
horizon). 

The predicted Rayleigh and Mie scattering flux of ZL, 
and respectively, along the line of sight to our 

targ et fi eld throughout our observations is shown in Fig¬ 
ures A6-^. In those Figures, we show the scattered light 
as a function of the above-the-atmosphere ZL flux in tar¬ 
get field at {a = 3.00h,i5 = —20.18d), 7(3h,-20d)(ZL). This 
removes the spectrum of the ZL and highlights the wave¬ 
length dependence of /scat ■ The predicted scattered flux is 
not symmetric about the zenith because the distribution of 
ZL over the sky is not symmetric: the scattered light will 
be smaller at the same airmass if the Sun is further below 
the horizon, i.e. the middle of the night. The scattered 
flux is therefore minimized near UT~4, when the field is 
still at low airmass and the brig htest regions of the ZL 
are below the horizon. In Figure AlO, we show the total 
combined effect of the atmosphere on the ZL flux received 
from the target field: 

/„et(A,t,x,ZL) = /gt,(A,t,x,ZL)±/g“,(A,<,x,ZL) 

-A3h,-20d)(A, ZL)(1 - e-"°-(^)^).(A6) 

Finally, from /net (A, ZL) we can calculate an effective ex¬ 
tinction for the ZL from our target field at the specific 
times at which our observations occurred. The effective 
extinction is defined by the equation 

/net(A,t,x,ZL) = (1 - e"»«(^’‘)x)/(3h,-20d)(A,ZL) (A7) 


The effective extinction is plotted in Figure All, and is 
specific to our target field, times of observation, observed 
extinction, geographic latitude and longitude, and alti¬ 
tude. 

The result which is applied to our ZL measurement from 
the modeling discussed here is Tgff (A, t), which corresponds 
to /net(A, t, Xj ZL) rather than /scat(A, t, Xi ZL). The virtue 
of this approach is that the abso lute flux accuracy of the 
adopted ZL over the sky (Fig. A4) does not affect our 
results; only the accuracy of the relative flux distribution 
over the sky matters. In the regions of the sky which 
dominate the scattering for our observations (solar elon¬ 
gations of e > 30°), the relative flux errors for the ZL over 
the visible hemisphere of the sky are <5% over large ar¬ 
eas (> 30°), and better on small scales. Such errors will 
propagate into final measurement of the ZL at the level of 
< 1% at high airmass, and < 0.4% at low airmass. Nev¬ 
ertheless, we note that the above-the-atmosphere value of 
the ZL from Levasseur-Regourd & Dumont (1980) agrees 
with our measurement in our target field to within 2%. 

To evaluate the accuracy of our calculated values of 
/net(ZL) (see Equation A7), we estimate that the uncer- 
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Wavelength (A) 


Fig. AlO.— The net effect of the atmosphere on the observed spec¬ 
trum of zodiacal light. 7net(ZL) is the flux received from the target 
field plus the scattered light coming from the entire visible hemi¬ 
sphere of the sky: Jnet(A, ZL) = 7scat(A, ZL) - 7(3;,^_20d)(ZL) (1 - 
g—7net(A, ZL) is plotted as a function of wavelength at, 
discrete times, as labeled. Flux units are the same as in Figures Aq 


Fig. All.— Each line shows the effective extinction for the zo¬ 
diacal light as a function of wavelength for our observations at the 
indicated UT. The effective extinction corresponds to the net loss of 
light relative to the above—the— aimpsphere flux of the zodiacal light 
in our held of view (see Figure AlO). For comparison with y^e total 
observed extinction derived from standard stars see Figure m. 


tainty in our scattering calculations is 8% at the low zenith 
angles (< 30°) where the bulk of our observations occur. 
This estimate is based on the comparisons between scat¬ 
tering models and atmospheric measurements presented in 
Green etal. (1971), Dave (1964), and Staude (1975), and 
is consistent with the uncertainties discussed in WvB67, 
Ashburn (1954), and Sekera & Ashburn (1953). The time- 
weighted average of /scat (ZL) over the course of our obser¬ 
vations is ~ 0.15 X /( 3 h.- 20 d)(^i ZL), so /scat(ZL) has an 
uncertainty of 1.2% of the ZL flux in our target field. The 
uncertainty in the observed extinction is much less than 
1% and adds negligibly to this error. See for further 
discussion of the accuracies of the zodiacal light measure¬ 
ment. 

We can also assess the accuracy of Teff(A) independently 
from our own data, as discussed in §5. Notice that /net(ZL) 
changes with time in a way which is only weakly dependent 
on wavelength. A consistent solution for the ZL with both 
wavelength and airmass will be strong conhrmation of the 
accuracy of the values for reff(A,t) calculated here. 

Integrated Starlight 

Unlike the scattered ZL, the scattering which results 
from integrated starlight (ISL) must be incorporated into 
our analysis of the observed night sky spectrum as an abso¬ 
lute flux value. We must therefore Hrst derive a spectrum 
for the ISL as a function of position over the sky. To do so, 
we have followed the method suggested by Mattila (1980a, 
1980b), which we briefly summarize here. 

The spatial and flux distribution of stars of all spectral 
types can be described by exponential distributions per¬ 
pendicular to the Galactic plane (in the z direction) and 
narrow Gaussian distributions in intrinsic magnitude. The 
mean emission per pc^ from stars of type i as a function 


of distance from the Galactic plane, z, can be written as 

Mz) = j,(0)e-l^l/'‘^ = A(0)10-°-4^^e-l^l/'‘S (A8) 

in which Di{0) is the number density of stars per cubic 
parsec in the plane, hi is the scale height of the verti¬ 
cal distribution, and Mi is the mean absolute magnitude 
of the spectral type i. The observed flux is also attenu¬ 
ated by interstellar extinction, which can be expressed by 
a two-component extinction law characterized by a total 
extinction ao(A) = ai(A) -I- 02 (A), with ai(A) : 02 (A) in 
the ratio 1.84 : 0.62 (Neckel 1965). The z-dependence of 
extinction can be written as 

a(z,A) =ai(A)/[l + (z/20)2]+a2(A)/[l + (z/100)"] (A9) 

for z given in parsecs (Neckel 1965, Neckel 1980). We 
find a good fit to the observed ISL by adopting standard 
values for ao(A) (Zombeck 1990), scaled to ao(U) = 1.5 
mag kpc“^. 

In cylindrical coordinates, the flux per unit solid angle 
(ergs s“^ cm“^ sr“^ A“^) from stars fainter than mg along 
the line of sight at Galactic latitude b can be expressed by 
the volume integral 

I^{b) = y 

^ 47rsm6 

(AlO) 

in which r is the distance along the line of sight from the 
observer in parsecs and fi is the spectral energy density of 
a star of type i in ergs s“^ cm“^ A“^. In the derivation 
of the above integral, the l/47rr^ loss of flux from each 
star along the line of sight has canceled with the in 

the volume integral, and we have changed variables using 
the relation r = z/sin 6. The lower limit of integration 
is simply the distance modulus for stars of each type cor¬ 
responding to the bright magnitude cut-off, ttiq, so that 
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5000 6000 

Wavelength (A) 


Fig. A12.— The mean integrated starlight (ISL) spectra as a 
function o f Oa. 1ar.tir latitude (as labeled) produced by the model 
described 5 4.3 


zq = in parsecs. Finally, the extinction 

from z to the observer is 

A^(z) = f a{z)dz (All) 

Jo 

= 20ai(A) arctan(z/20) + 100a2(A) arctan(z/100). 


Using 33 individual stellar types described by the param¬ 
eters Mi, Di, and hi from Wainscoat etal. (1992), we ob¬ 
tain integrated spectra which agree with the observed star 
counts at V and B (Roach & Megill 1961, see also Allen 
1973) to Too = 6 U mag at |5| > 5° to better than 10%, 
which is more than adequate for our purposes. The spec¬ 
tral energy densities for each stellar type, /i(A) were ob¬ 
tained from the STScI archive (Jacoby, Hunter, & Chris¬ 
tian 1984) and have a resolution of roughly dA. We include 
stars by type with mg < 6 U mag from the SAG star cat¬ 
alog by hand. We felt this was necessary as the statistical 
variation in stellar density on small scales around the solar 
neighborhood can have a significant impact on the accu¬ 
racy of the model, while variations are apparently averaged 
out in stellar populations at large distances. In Figure 


A12, we plot the total integrated starlight (ISL) with no 


bright magnitude cut off at 0° < |6| < 90°. The total flux 
at |6| = 90° is roughly 20x10“® ergs s“^ cm“^ sr“^ A“^, 
while the flux near the plane is as high as 300x10“® ergs 
s“^ cm“^ sr“^ A“^. Interstellar extinction limits the ISL 
flux in the plane in our model, probably more than is ap¬ 
propriate. However the flux rises rapidly even 1° degree 
above the plane to more realistic values. The limited sky 
area at 5 = 0° precludes this from impacting the accuracy 
of the mo dels. 

Figure A15 shows the total scattered ISL flux due to 
Rayleigh and Mie scattering which contributes to observa¬ 
tions of our target field at the beginning of the observing 
night (UT=2.0 hr). The total flux (~ 12 x 10“® ergs s“^ 
cm“^ sr“^ A“^) is roughly 12% of the ZL flux above the 
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Fig. A13.— Integrated starlight over the sky from star counts at 
V. Flux units are 1x10“^ ergs s“^ cm“^ sr“^ A~^. 



4000 4200 4400 4600 4000 5000 

Wavelength (A) 


Fig. A15.— The upper plot shows the predicted spectrum of the 
scattered ISL along the line of sight to the EBL field at UT=2.0 hr 
for our observations. Units are ergs s~^ cm“^ sr“^ A~^. The solar 
spectrum is also shown, scaled to the same flux and offset, to allow 
visual comparison of the spectral features. The lower plot shows the 
ratio of that spectrum to the solar spectrum, normalized at 4600A, 
the center of our observed wavelength range. Note that the total 
flux from the scattered ISL in this case is < 5% of the total flux of 
the combined zodiacal light and airglow, and the spectral features 
are weaker in the scattered ISL. 


atmosphere in our target field. The ISL flux in the ob¬ 
served sky spectrum will impact our measurement of the 
ZL flux only if the ISL spectrum has the same features as 
the sol ar (z odiacal light) spectrum. In the lower plot of 
Figure A15 , we therefore plot ratio of the /scat(A,t,x,ISL) 
to the solar spectrum, normalized at 4600A. It is clear 
from this plot that the scattered ISL and solar spectra dif¬ 
fer by 5-10% at > 4500A, but by a factor of 3-5 in the 
stren gth of spectral features at less than 4500A. In Figure 
AI6|, we plot the ratio of /scat (A, t, y, ISL) throughout the 
night to /scat(A,t,x,ISL) at UT=Ao hr. From this plot it 
is clear that strength of spectral features changes only very 
weakly throughout the night, by < 1% over the majority 
of the spectrum and by < 4% at 3900-4000A (Cal H & 
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UT= 2.0 


UT= 3.1 




UT= 7.5 



Fig. A14.— Integration pattern for the calculation of scattered ISL flux. Each plot is an Aitoff projection of the sky in Galactic coordinates. 
The center of the Galaxy is at the center of each plot (^ = 0°, 6 = 0°). The Galactic plane within 30 degrees of the Galactic center is marked 
by the thick line. Our target field is indicated by the square. The Galactic coordinates of the zenith can be seen as the “bullseye” center of 
the integration pattern. The coordinates of the local horizon are shown by the edge of the integration pattern. The Galactic plane is running 
along the horizon at the start of the night (TT T=9 0 hr), and is perpendicular to the horizon at the end of the night (UT=7.5 hr). The Galactic 
center is never above the horizon. See Figure A5 for further discussion. 



Wavelength (A) 

Fig. A16.— The ratio of the predicted scattered ISL spectrum 
at the indicated UT on the night of our observations, relative to 
the predicted spectrum at UT=2.0 hr. The plot demonstrates that 
relatively small changes (< 4%) occur in the strength of individual 
spectral features over the night, while changes in broad-band color 
and mean flux are significant. 


K). The consistency of our ZL measurement (§^) over the 
full wavelength range 3900-5100A is, therefore, a strong 
test of the accuracy of the predicted contribution of scat¬ 
tered ISL. As discussed in the predicted scattered ISL 
flux is entirely consistent with our observations. See for 
further discussion. 


Obviously, the model we describe above makes no al¬ 
lowance for variation in the ISL with Gala ctic longitude. 
For comparison, we show in Figure A13 an Aitoff pro¬ 
jection of the ISL from star counts over the sky, which 
shows that the ISL has only minor dependence on longi¬ 
tude at 5 > 20°. At lower latitudes where the variation 
is greater with longitude, the spectroscopic model we em¬ 
ploy does give a good approximation to the average ISL. 
Because the contribution to the scatterin g com es f rom a 
wide spread in longitude (compare Figures A14 and AI3), 
the mean value is adequate for our purposes. To test this, 
we ran simulations in which we maintained the mean ISL 
flux with latitude, but varied the ISL flux with longitude 
by ±25%. We find that the total scattered ISL is affected 
by less than 9% throughout the night due to longitudinal 
variations around the mean. 

The mean flux in our models for the ISL is consistent 
with the star counts of Roach & Megill (1961) to within 
±10% at both V and B. As in the previous section, we 
estimate that the uncertainty in our scattering calculations 
is 8% at the low zenith angles (< 30°) where the bulk of 
our observations occur. Combining these, we estimate the 
uncertainty in /scat(A, t, y, ISL) to be 13%. As the relative 
mean strength of spectral features in starlight is 0.6-4% of 
the total ZL flux in our target field at the 4000 — 5200A, 
this corresponds to an uncertainty of < 0.5%. 

Any significant errors in our model, either in mean flux 
as might be caused by longitudinal variations in the ISL, 
or in the spectral energy distribution, would show up as in¬ 
consistencies in the solution for the ZL flux as a function of 
wavelength. Furthermo re, er rors would be worst at higher 
airmass, where Figure A14 shows that the low-galactic- 
latitude sky has a greater impact on the scattered ISL, 
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the mean flux is greater, and the stellar-type mix is more 
sensitive. No such variations with wavelength are found, 
as we have discussed in §|[ See ^ for further discussion 
of the accuracies of the zodiacal light measurement. 
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